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Abstract. The original continuous-time "goldfish" dynamical system is characterized by 

■ two neat formulas, the first of which provides the N Newtonian equations of motion of this 
| dynamical system, while the second provides the solution of the corresponding initial-value 

problem. Several other, more general, solvable dynamical systems "of goldfish type" have 
^ \ been identified over time, featuring, in the right-hand ("forces") side of their Newtonian 

■ equations of motion, in addition to other contributions, a velocity-dependent term such as 
that appearing in the right-hand side of the first formula mentioned above. The solvable 
character of these models allows detailed analyses of their behavior, which in some cases is 
quite remarkable (for instance isochronous or asymptotically isochronous) . In this paper we 
introduce and discuss various discrete-time dynamical systems, which are as well solvable, 
which also display interesting behaviors (including isochrony and asymptotic isochrony) and 
which reduce to dynamical systems of goldfish type in the limit when the discrete-time 
independent variable I = 0,1,2,... becomes the standard continuous-time independent 
variable t, < t < oo. 
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q . 1 Introduction 

i— i ' 

The original "goldfish" dynamical system [IJ [2] is characterized by the system of N Newtonian 
equations of motion 

N 2z n z k 

z n = ^2 — ^ — , n = l,...,N, (1.1a) 

and by the following neat prescription yielding the solution of the corresponding initial-value 
problem: the N values of the dependent variables z n = z n (t) at time t are the iV solutions of 
the algebraic equation (for the unknown z) 

£-«°Li (i.ib) 

^z-z fc (0) t 

i.e. the N roots of the polynomial equation of degree iV in the variable z that obtains by 

N 

multiplying this formula by the polynomial Yl [ z ~ z j(ty]- 

3=1 



*This paper is a contribution to the Proceedings of the Conference "Symmetries and Integrability of 
Difference Equations (SIDE-9)" (June 14-18, 2010, Varna, Bulgaria). The full collection is available at 
|http://www.emis.de/journals/SIGMA/SIDE-9.html| 
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Notation 1.1. Here and hereafter N is an arbitrary positive integer (generally N > 2), su- 
perimposed dots denote differentiations with respect to the independent variable t ("continuous 
time"), and the N dependent variables z n = z n (t) may be interpreted as the coordinates of N 
point-like unit-mass particles - hence z n denote their velocities and z n their accelerations, consis- 
tently with the interpretation of (jl . la|) as a set of Newtonian equations of motion with velocity- 
dependent forces. The indices - such as n, m, j, k - generally run from 1 to N (below, as 
a convenient reminder, we often indicate this explicitly; as well as the exceptions to this rule). 

Hereafter we denote as "dynamical system of goldfish type" any dynamical system charac- 
terized by Newtonian equations of motion featuring in their right-hand sides - which have, in 
the Newtonian context, the significance of "forces" - a velocity-dependent term such as that 
appearing in the right-hand side of (|l.la|) (of course in addition to other terms). Let us also 
emphasize that the dynamical system (jl.ip is the simplest model belonging to the Ruijsenaars- 
Schneider integrable class [313]. 

For instance a simple extension of the above model (reducing to it for uj = 0) is characterized 
by the Newtonian equations of motion: 

z n = (1 — 2a)iujz n + a(a — \)uj 2 z n 

^ 2(z n + iaujz n )(z m + iaujZm) 
+ 2^ ~ > n = l,...,N. (1.2a) 

-. , Zn Z m 

The corresponding solution of the initial-value problem is again given by a simple rule: the N 
values of the dependent variables z n = z n (t) at time t are related by the formula 

z n {t) = ( n (t)exp(-iau;t), n = l,...,N, (1.2b) 

to the N solutions Cn{t) of the algebraic equation (for the unknown £) 



v 



k=l 



E ffc(O) + iaujz k (0) _ iu 
C - 2 fc (0) ~ exp(iwt) - 1 



(which again gets transformed into a polynomial equation of degree N in £ after multiplication by 

N 

the polynomial Y\ [C ~ z j(0)]); or, equivalently but more directly, the N values of the dependent 
i=i 

variables z n = z n (t) at time t are the N solutions z n (t) of the algebraic equation (for the 
unknown z) 

ifc(O) + iauzkjO) _ iu exp(iaujt) 
z — Zfc(0) exp(— iauit) exp(iut) — 1 

Notation 1.2. Here and hereafter i is the imaginary unit, i 2 = —1, to is an arbitrary constant 
(dimensionally, an inverse time), and a is an arbitrary (dimensionless) constant. Clearly - unless 
both u and aoj are both imaginary, Re(cj) = Re(acj) = - the time-evolution of this system 
takes place in the complex z-plane, i.e. the dependent variables z n = z n {t) are complex; but 
it may as well be viewed as describing the evolution of N point-like particles moving in the 
real xy-plane - whose positions at time t are characterized by the (real) Cartesian coordinates 
x n = x n 

(t), Vn = Vn{t) ~ by setting z n (t) = x n (t) + iy n (t); and one of the remarkable features 
of the resulting real model is the possibility to write its Newtonian equations of motion in 
covariant - i.e., rotation-invariant - form, see Chapter 4 of [3]. Hereafter we generally refer to 
this model, and its generalizations, see below, in their complex versions. 
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Remark 1.1. Clearly for to = this model, (jl.2p . reduces to the previous model (jl.ip . For to ^ 
the time evolution of this model, (|1.2|) . depends mainly on the values of the two constants a; 
and ato, as displayed by its solution, see (|1.2d[) . If both these constants are real, Im(w) = 
Im(aw) = (hence as well Im(a) = 0), the time evolution of this model is confined, indeed 
completely periodic if the real number a is rational, while if a is irrational it is multiply periodic, 
being a nonlinear superposition of two periodic evolutions with the two noncongruent periods 
T = 2tv/ \lo\ and T/a. Note that these outcomes obtain for generic initial data: hence, for a 
rational, a = q/p with q and p coprime integers (and p > 0), this system is isochronous, its 
generic solutions being completely periodic with period pT - or possibly with a period which 
is a, generally small, integer multiple of pT: indeed, when the equation (|1.2dp is itself periodic 
with period pT, the unordered set of its N roots is clearly periodic with the same period pT, 
but the periodicity of the time-evolution of each individual coordinate z n (t) may then be a, 
generally small, integer multiple of pT due to the possibility that different roots get exchanged 
through the time evolution (for a discussion of this phenomenology - including a justification 
of the assertion that the relevant integer multiple of pT is generally small - see [5]). 

On the other hand, if u> is real but a is imaginary, say ato = i^ with 7 real and nonvanishing, 
then clearly in the remote future - i.e., as t — > 00, and up to relative corrections of order 
exp(— |7|i) - all the N coordinates z n {t) tend to the origin, z n {oo) = 0, if 7 < 0, while if 7 > 
they all diverge (see (11.2bj) and (jl.2c|) ). 

If instead to is not real, Im(w) 7^ 0, then in the remote future (i.e., as t — > 00, and up to 
relative corrections of order exp(— | lm(cj)\t)) the N solutions of (jl.2c() become asymptotically, if 
Im(u;) > 0, the N solutions £ n = Cn(oo) of the time-independent polynomial equation of order N 
in £ 

^ ifc(O) + iaujz k (0) _ _. 



k=l 



C - **(0) 



while if instead Im(oj) < the equation (|1.2c[) becomes, in the remote future, the time- 
independent polynomial equation 



ifc(O) + iaujz k (0) 

k= 



0. 



hence N — 1 of the solutions of (jl.2cp tend asymptotically to the N —1 solutions of this equation 
(polynomial of degree N — 1 in £) and one of them approaches asymptotically the diverging 
coordinate 



N 



Casy(i) = exp(iwi) y^[ifc(0) + iauz k (0)]. 



k=l 

Note that this implies (see (|1.2bj) ) that, if Im(a;) > but aw is real, ato = p with p real and 
nonvanishing, then the model f)1.2|) is asymptotically isochronous, its generic solutions beco- 
ming, in the remote future, completely periodic with period 2n/\p\, up to corrections vanishing 
exponentially as t — > 00 (for a more detailed discussion of the notion of asymptotic isochrony 
see Chapter 6, entitled "Asymptotically isochronous systems", of [6]). 

For to = (i.e., when the model (ll.2p reduces to (II. ip ) it is possible to restrict consideration 
to real dependent variables z n , but even then it is more interesting not to do so, so that the 
time evolution takes place in the plane rather than on the real line: see the remarkable behavior 
of this dynamical system in this case ("the game of musical chairs"), as detailed in Section 4.2.4 
of [1]. Hence let us reiterate that we always consider the dependent variables z n to be complex 
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numbers, both in the continuous-time case, z n = z n (t), < t < oo, and (see below) in the 
discrete-time case, z n = z n (£), £ = 0, 1, 2, ... . 

Another large class of solvable dynamical systems "of goldfish type" is characterized by the 
equations of motion 

N 

z n = aiz n + a 2 + a 3 z n - 2(N - \)a^z 2 n + ^ (z n - z m )~ l [2z n z m 

m=l, mj^n 

+ (o 5 + a 6 z n )(z n + z m ) + a 7 z n (z n z m + z m z n ) + 2(a 8 + a 9 z n + a w zl + a 4 ^)] , 
n = l,...,N, (1.3) 

featuring 10 arbitrary constants (see [U equation (2.3.3-2)]). In this case the solvability is 
achieved by identifying the ./V dependent variables z n (t) with the N roots of a time-dependent 
polynomial ip(z, t) of degree N in z satisfying a linear second-order PDE in the two independent 
variables z and t. 

For an explanation of the origin of the name "goldfish" attributed to these models see Sec- 
tion l.N of [B] and the literature cited there. In this book [B] (see in particular its Section 4.2.2, 
entitled "Goldf ishing" , and the papers referred to there) several other solvable models "of 
goldfish type" are reported, including isochronous ones (i.e., models featuring solutions which 
are completely periodic with a period independent of the initial data). A few additional models 
of goldfish type have been identified more recently [HE]. 

The most remarkable aspect of these dynamical systems is their solvability, namely the possi- 
bility to solve their initial- value problems by algebraic operations, amounting generally to finding 
the N eigenvalues of an N x N explicitly known time-dependent matrix (see below), or equiva- 
lently to finding the N roots of an explicitly known time-dependent polynomial of degree N (see 
for instance (jl.lbp and f)l .2d[) ) . Quite interesting is also the identification of multiply periodic, 
completely periodic, or even isochronous or asymptotically isochronous cases. 

In the present paper we present various discrete-time dynamical systems "of goldfish type" , 
so denoted because all these models reduce, in the limit when the discrete-time independent 
variable £ = 0,1,2,... becomes continuous, to continuous-time dynamical systems of goldfish 
type. All these models are moreover solvable, i.e. the solution of their initial value problems can 
be achieved by finding the N eigenvalues z n {£) of N x N matrices explicitly known in terms 
of the initial data and of the discrete-time independent variable £; or equivalently by finding 
the N roots z n {£) of a polynomial, of degree N in the complex variable z, as well explicitly 
known in terms of the initial data and of the discrete-time independent variable £. Some of 
these models feature interesting behaviors, even isochrony or asymptotic isochrony. Two of 
these models (see Subsection 12.11 and I2.2|) were treated in the paper [10] . which has not been 
published because - after it was submitted for publication but before getting any feedback - 
new solvable models were identified and it was therefore considered preferable to report all these 
models in a single paper, this one. The main properties of each of these discrete-time models 
are reported in Section El and proven in Section [3j These properties include the display of the 
equations of motion of these discrete-time models, the solution of their initial-value problems, 
a terse discussion (for the first three models) of their behavior including the possibility that 
for special values of some of their parameters they possesses periodic or multiperiodic solutions 
or even display isochrony or asymptotic isochrony, and some mention of their continuous-time 
limits. Section U] entitled "Outlook" concludes the paper: in it a general framework is outlined 
which might allow the identification of additional solvable discrete-time models. And some 
mathematical developments are confined to two appendices. 

These findings are congruent with the recent surge of interest for discrete-time evolutions - 
in particular, such evolutions which are in some sense integrable or even solvable. Given the 
large body of research devoted to these topics over the last two decades, our reference to the 
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relevant literature shall be limited to citing the following surveys: [11[ [T2| [T3l [HI |15| , But 
a special mention must be made of the seminal papers by Nijhoff, Ragnisco, Kuznetsov and 
Pang, see [16] as well as the earlier paper [17], where discrete-time versions were introduced of 
the well-known integrable "Ruijsenaars-Schneider" and "Calogero-Moser" dynamical systems; 
as well as of the paper by Suris [18], which treats specifically a discrete-time version of the 
original goldfish model. Some results of these papers refer to models whose equations of motion 
feature trigonometric/hyperbolic or even elliptic functions, and are therefore more general than 
those treated in the present paper, whose equations of motion only feature rational functions 
(see below); on the other hand the findings reported below include more general models than 
those previously treated, demonstrate the solvability of these models by an approach somewhat 
different from those previously employed, and, most significantly, display the possible emer- 
gence of remarkable phenomenologies - including periodicity and even isochrony or asymptotic 
isochrony, see below - not previously identified for this kind of discrete-time dynamical systems. 

Let us also mention that the approach developed below also allows to identify and investi- 
gate discrete-time variants of another class of solvable continuous-time dynamical systems, the 
prototype of which is characterized by the Newtonian equations of motion 



(with c an arbitrary constant), instead of (jl.lap . But in this paper we merely indicate, at the 
appropriate point, how to proceed in this direction, postponing a complete treatment of this 
development to a separate paper. 

And let us finally pay tribute to Olshanetsky and Perelomov who were the first to show, more 
than 35 years ago, that the time-evolution of a nontrivial many-body system could be usefully 
identified with the evolution of the eigenvalues of a matrix itself evolving in a much simpler, 
explicitly solvable, manner: see [19] and their other papers referred to in Section 2.1.3.2 of [4], 
entitled "The technique of solution of Olshanetsky and Perelomov" . The present paper extends 
their approach to the discrete-time context. 

2 Results 

In this section we report the main results of this paper; they are then proven in Section [3l 

Notation 2.1. Hereafter the dependent variables are indicated again as z n , but they are now 
functions, z n = z n (£), of the discrete-time variable I taking the integer values I = 0, 1, 2, ... ; and 
superimposed tildes indicate generally a unit increase of the independent variable £, for instance 
Zn = z n (£ + 1), z n = z n {£ + 2). Hereafter 

&nm is the standard Kronecker symbol, 5 nm — 1 if 
n = m, 5 nm = if n ^ m, and underlined quantities are iV-vectors, for instance z=(z±,..., zn). 
For the remaining notation we refer to Notation II. 1[ see above, and to specific indications given 
case-by-case below. 

As reported in this section and explained in Sections [3] and [U the solvable models considered 
in this paper generally feature three equivalent versions of the second-order "equations of mo- 
tion" characterizing their evolution in discrete-time. The treatment of the first model given in 
Subsections 12.11 and 13.11 is somewhat more detailed than that provided for the other models in 
the subsequent subsections where, to avoid repetitions, we often refer to the treatments provided 
in Subsections 12.11 and 13.11 And already in this section, as well as in Section [3] we often take 
advantage - to simplify the presentation of some results - of identities and lemmata collected 
in Appendix [Al 
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2.1 First model 

The first model is defined by the following second-order discrete-time equations of motion: the 
N values of the twice- updated variables z n = z n (£ + 2) are given, in terms of the 2N values of 
the variables z m = z m (£), z m = z m (£ + l), by the N roots of the following (single) algebraic 
equation in the unknown z, 



N 

E 

k=l 



z k - az k 



z 



az k 



N 

n 



Z k Zj 



1. 



(2.1a) 



which clearly amounts to a polynomial equation of degree N in this variable z (as it is imme- 

N 

diately seen by multiplying this equation by the polynomial FJ (z — az m )). Here and below a is 

m=l 

an arbitrary (dimensionless, nonvanishing) constant. A neater version of this formula is easily 
obtained by multiplying it by a and by then using the identity ()A.10|) with rj n = az n , Q n 



n 



1, . . . , N. It reads 



N 

n 



z — a Zj 



z — aZj 



1 + a. 



a Zr, 



(2.1b) 



An equivalent formulation of this model is provided by the following system of N polynomial 
equations of degree N for the twice-updated coordinates z n = z n (£ + 2): 



JV 

E 

k=l 




N-l 



n 



,N. 



(2.2a) 



Again, a neater version of this formula is easily obtained by dividing it by a and by then using 

1, . . . , iV. It reads 



the identity ()A.10|) . now with z = a 2 z n , rj n = a z n , ( n = z n , n 




(1 + a)a 



N-l 



n 



N. 



(2.2b) 



And a third, equivalent version of this model is provided by the following system of N 

+ 2): 



polynomial equations of degree N for the twice-updated coordinates z n = z 

N 




N-l 



n 



1,...,N. 



(2.3) 



The similarities and differences among these three sets of "equations of motion", (|2.ip . (I2.2h 
and (|2.3p . are remarkable: let us reemphasize that they in fact yield the same evolution in 
discrete-time of the N coordinates z n = z n (£). Particularly remarkable is their similarity in the 
special a = 1 case, when the 3 versions (|2.1bp . (|2.2bp and (j2.3|) of the equations of motion read 
as follows: 




n 



n 



Discrete-Time Goldfishing 



7 




1. 



n 



1,...,N. 



The last of these three systems coincides with equation (1.8) of [18J . 

Remark 2.1. This model, see (|2.ip . (|2.2p and f|2.3|) - as the original goldfish model, see (|1.2p - 
is invariant under an arbitrary rescaling of the dependent variables, z n => cz n with c an arbitrary 
constant; including the special case c = exp(i7) with 7 an arbitrary real constant, corresponding 
to an overall rotation around the origin in the complex z-plane. 

The solution of the initial-value problem for this model is given by the following 

Proposition 2.1. The N values z n (£) of the dependent variables at the discrete time i are the 
N eigenvalues of the N x N matrix 



U nm (£) = 5 nm z n (0)a e + v m (0)- 



1 



n, m 



,N 



(2.4a) 



with 



V m = V m (z,i) 



N 

3=1 



az Tl 



N 



]1 [a(zj - z m )] 

i=lj j+rn 



m 



(2.4b) 



where of course v m (0) indicates the value ofv m (z,z) corresponding to the initial data z = z(0), 
5 = 5(0) =l(l). 

A neater, equivalent formulation of this finding - obtained from (|2.4p via Lemma IA.4I with 
Cn = z n (0)a e and r\ m = v m (0)(a — l)/(a — 1) - states that the N coordinates z n {€) are the N 
solutions of the following algebraic equation in z: 




z k (l) - az k (0) 



z - a e z k {0) 



N 

n 

J'=l, 3¥=k 



Zj{\) - az k (0) 
azj(0) — az k (0) 



1 



And another, even neater, equivalent formulation - obtained from this via the identity (jA.lOp 
with z replaced by za 1 ~ i , r]k = az k (0), Q = %(0) = Zj(l) - states that the coordinates z n {t) are 
the N solutions of the following algebraic equation in z: 



N 

n 

k=l 



l z k (l) 



a £ z k (0) 



J-i 



1 



(2.5) 



The last two equations become of course polynomial equations of degree N in z after multiplica- 

N 

tion by the product Y\[ z ~ a z j(®)]- 
3=1 

These formulas are also valid for a = 1 (by taking the obvious limit, i.e. replacing (a? — 
l)/(a — 1) with p). If instead \a\ < 1, then clearly for all (positive) values of I the matrix U{€) 
is bounded and U nm (oo) = v m (0)/(l — a); hence for all values of i the iV coordinates z n {€) are 
bounded and N — 1 of them vanish as i — > 00 while one of them tends to the value 



iV 



-asy 



( l_ a )-i^ fc(0 ) 



1 



k=l 



N 
k=l 



z k (l) - azfc(0)], 
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see (pnbj) and the identity (|ATT|) (with n k = a z k (0), (j = %(0) = Zj(l)). If a ^ 1 but it has 
unit modulus, 

a = exp(27riA) (2-6) 

with A real and not integer, then clearly the matrix U{£) is again, for all values of £, bounded; 
and if moreover A is a (strictly, i.e. non integer) rational number, 

A = — , a = exp( — —J, (2.7) 

with K and L two coprime integers and L > 1, then clearly the matrix ?7(£) is periodic with 
period L, 

f/(/ + L) = U(l), 

hence the (unordered) set of its iV eigenvalues z n (£) is as well periodic with period L. This 
shows that in this case, see (|2.7p . the discrete-time goldfish model, see (|2.ip or (| 2 . 2 j) or (|2.3p . 
is isochronous. On the other hand if A is real and irrational, then clearly the time evolution of 
this discrete-time dynamical system is not periodic: indeed, while the right-hand side of (|2.4a|) 
(with (12. 6p and A real and irrational) is periodic (with unit period) as a function of the real 
variable r = X£, clearly it is not periodic as a function of the variable £ taking the integer 
values t = 0, 1, 2, . . . . (We made this analysis, for convenience, referring to the coordinates z n (£) 
as the eigenvalues of U(£), see (12. 4ft : of course an analogous discussion could be made on the 
basis of the alternative identification of the coordinates z n (£) as the N roots of the polynomial 
equation (12. 5j) - whose similarity with (|1.2dj) is in any case to be noted, see below.) 

To explore the transition from the discrete-time independent variable I to the continuous-time 
variable t one makes the formal replacements 

I => -, £+1 — , £ + 2 => t -^, (2.8a) 

e e e 

a 1 - iwe, (2.8b) 
and (with a slight abuse of notation) 

ZnW ^ > Zn(*), 

z n (£) =z n (£ + l) =► z n (t) + ez n (t) + —z n {t) + 0(e 3 ), 

i n {£) = Zn {£ + 2) z n (t) + 2ei n (t) + 2£ 2 z n (t) + 0(e 3 ), (2.8c) 

with e infinitesimal. It is then a matter of standard, if a bit cumbersome, algebra, to verify 
that the insertion of this ansatz, see ()2.8bp and (|2,8cp . in (j2. lb|) or (|2.2b|) or (|2.3p yields a trivial 
identity to order e = 1, while to order e it reproduces (|1.2aj) with a = 1, reading 

A 2(i n + iwz n )(i m + iwz m ) 
z n = -zwz n + 2^ ' n = l,...,N. (2.9a) 

m=l, m^tn 

Likewise, the discrete-time solution formula (|2.5[) becomes, in the continuous-time limit, 

ifc(O) + iuz k (0) _ iuj _ wexp(wt) gb ^ 

^-^ z — Zfc(0) exp(— iwt) 1 — exp(— iut) exp(iwt) — 1 ' 

which coincides with (|1.2dp with a = 1. A terse outline of the derivation of these results 
is provided at the end of Subsection 13.11 To higher order in e one would obtain additional 
relations satisfied by the solution z n (t) of this continuous-time goldfish model, which might 
alternatively be obtained by differentiating its equations of motion (|2.9al) . 
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Remark 2.2. Clearly, for ui real and nonvanishing, this continuous-time model, (|2.9ap . is 
isochronous: see ()2.9bp and/or Remark II .11 This is consistent with the fact that the limiting 
replacement ()2.8b|) can be considered to obtain from (|2.7p - entailing isochrony of the discrete- 
time model - by identifying eoj with 2irK/L in the context of the replacement (see (|2.8p ) of 
the unit interval in the discrete-time model with the infinitesimal time interval e to make the 
transition to the continuous-time case. 

Remark 2.3. At every step of the discrete-time evolution the N values of the twice-updated 
variables z n = z n (£ + 2) are given, in terms of the N unupdated variables z m = z m (£) and the N 
once-updated variables z m = z m (£ + 1), as the N roots of a polynomial, of degree N in its argu- 
ment z, whose coefficients are explicitly defined in terms of the 2N unupdated and once-updated 
variables: see (|2.1|) or (|2.2|) or (|2.3|) (hereafter - within this important Remark l2.3l - we generally 
identify, for simplicity, this model only via the version (|2.1|) of its equations of motion). Hence 
at every step of this discrete-time evolution the unordered set of N twice-updated variables z n 
is uniquely determined, but not the value of each of them. This implies a qualitative difference 
among the continuous-time respectively the discrete-time evolutions described by the equations 
of motions (|2.9aj) (or, more generally, (|1.2aj) and (J1.3P ) respectively by (|2.ip : in contrast to 
the continuous-time case, the discrete-time evolution (12. ip is only deterministic in terms of the 
unordered set of N coordinates z m (£), but not for each individual coordinate z n {£). Indeed the 
continuous-time Newtonian equations of motion, see for instance (|1.2ap . determine uniquely the 
value of the acceleration z n (t) of the n-th moving point in terms of the iV positions z m (t) and 
the N speeds z m (t) of all moving points; and correspondingly, while the solution formula (jl.2d|) 
determines only the unordered set of N values z n (t) as the N roots of a polynomial of degree N, 
the value of each individual coordinate z n (t) gets then uniquely determined by continuity in 
the time variable t. This latter mechanism to identify uniquely the value of the coordinate of 
each moving point is instead missing for the discrete-time evolution (|2.ip . On the other hand 
it is clear that there are appropriate ranges of values of the parameter a and of the 2N initial 
data z m (0), z m (0) = z m (l) - with a sufficiently close to unity, the ./V initial coordinates z m (0) all 
sufficiently well separated among themselves, and each z m (0) = z m (l) sufficiently close to the 
corresponding z m (0), see (|2.8p - which cause the evolution yielded by the discrete-time goldfish 
model (12. ip to mimic closely that yielded by the continuous-time goldfish model (12.9p . provided 
at every step of the discrete-time evolution the appropriate identification is made of the value 
of each twice- updated coordinate z n = z n {£ + 2) (among the unordered set of N values yielded 
by the discrete-time equations of motion) by an argument of contiguity with z n = z n (£+ 1) and 
z n = z n{£)] an d likewise an appropriate identification is made by contiguity of each coordinate 
z n (£ + 1) with the corresponding coordinate z n {£) (among the unordered set of N values yielded 
by Proposition ^. ip - these arguments of contiguity taking the place of the continuity of z n (t) as 
function of t applicable in the continuous-time case. But the contiguity argument breaks down if 
the positions at time £ of two different points, z n (£) and z m (£) with n ^ m, get too close to each 
other, corresponding to a quasi-collision, or even coincide, corresponding to an actual collision; 
which is however not featured by the generic solution of the discrete-time model f)2. 1 1) - nor of 
the standard goldfish models (jl.2ap or (jl.3p - clearly emerging only for a set of initial conditions 
z n (0), z(0) = z n (l) having unit codimension in the 2A r -dimensional (complex) phase space z, z. 

This important remark is applicable to all the discrete-time models considered below, al- 
though it will not be repeated. 

Remark 2.4. Several of the formulas written above (in this section) simplify somewhat via the 
following replacement of the dependent variables: 



z n {£) => a e z n {£), n = l,...,N. 
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In particular the 3 equivalent versions (|2.1b|) , (|2.2bp and (|2.3p of the discrete time equations of 
motion are thereby reformulated to read 

JV 

n 

N 

n 

iV 

n 

i=i 



~ n 


~ z i\ 


= l + a 


Z n 


~h) 




Z i 


- Z n \ 


_ 1 + a 


h 


- z n j 


a 










~ Zn j 





n = l,...,N, 
n = l,...,N, 



and correspondingly the formula (|2.5p providing the solution of the initial-value problem reads 



N 

n 

k=l 



z n (£) - a 1 z k {l) 
z n (l) - z k (0) 



a e - 1 



n = l,...,N. 



Somewhat analogous remarks are applicable to all the discrete-time models considered below; 
their explicit implementation is left to the interested reader. 



2.2 Second model 

In this subsection we treat rather tersely a discrete-time dynamical system that generalizes 
the discrete-time goldfish model described in the preceding Subsection 12.11 This generalization 
amounts to the presence of an additional free parameter, b: indeed, for b = one reobtains the 
model treated in the preceding Subsection 12.11 (hence in this subsection we assume that b does 
not vanish, 67^0). 

The three equivalent versions of the equations of motion of this model read as follows. The 
first version identifies the twice updated coordinates z n (£-\-2) as the N solution of the following 
equation in z (amounting to the identification of the N roots of a polynomial of degree N in 
this variable): 

y[(~z k -az k \ fl_±bh\ yr W~z k - azA fl + 6* 3 -/«\1 I = 
The second and third versions consist of the following two systems: 

n = l,...,N; (2.10b) 
n = l,...,N. (2.10c) 

Remark 2.5. This model - as the original goldfish model (|1.2aj) . and as the model treated 
above, see Remark 12.11 - is invariant under a rescaling of the dependent variables, z n cz n 
with c an arbitrary constant; but only provided the parameter b is also rescaled, b =>■ b/c. 

The solution of this model is provided by an analog of (the first part of) Proposition I2.H 
reading as follows: 




.JV-l 



N-l 



'1 + bz n 
(1 + bz n 
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Proposition 2.2. The N values z n (£) of the dependent variables at the discrete time £ are the 
N eigenvalues of the N x N matrix 

U(£) = U(0)[al + bV(0)} e + V(0)[(a-l)I + bV(0)}~ 1 {[aI + bV(0)] e - I}, (2.11a) 
where again (see (|2.4D ) 

U(0) = diag[z n (0)], U nm (0) = <W„(0), (2.11b) 
while the N x N matrix V(0) is now defined componentwise as follows: 

[V(0)]nm = : "T^L . n,m = l,..., N, (2.11c) 
1 + bz m (0) 

with the quantities v m (0) defined again as in Subsection \2.1\ (see ()2.4bp and the sentence following 
this formula). Let us recall that I is the N x N unit matrix (whose presence in (|2.11a|) . however, 
might well be considered pleonastic) . 

As evidenced by a comparison of (|2.11ap with (|2.4ap . the behavior of the solutions of this 
model (with b ^ 0) are less simple than those of the model discussed in the preceding Subsec- 
tion [2TT1 In particular a confined behavior emerges only, see (|2.11ap . from initial data z n (0), 
z n (0) = z n (l) implying, via (|2.11cp and (|2.4bp . that all the N eigenvalues of the N x N matrix 
al + bV(0) have modulus not larger than unity, reading exp(— q n + 2-7rir n ) with the numbers q n 
and r n real and the N numbers q n nonnegative, q n > 0. If moreover the N numbers q n all 
vanish and the N numbers r n are all rational, the behavior is periodic (but not isochronous, 
since these numbers, q n and r n , generally depend on the initial data; see (|2.11c|) and (|2.4bp ). 
While, if some of (but not all) the N numbers q n are positive, and none is negative, then the 
phenomenology we just described (corresponding to the g n 's all vanishing) emerges only asymp- 
totically, as £ — ^ 00, up to corrections of order exp(— q£) with q the smallest of the nonvanishing 
numbers q n - provided all those r n 's are rational whose corresponding q n vanish. 

Let us finally mention that, also for this second model, a transition from the discrete-time 
independent variable £ to the continuous-time variable t can be performed (as tersely outlined 
at the end of Subsection I3.2[) : but the continuous-time goldfish- type model obtained in this 
manner turned out to be, to the best of our knowledge, new, hence it seemed appropriate to 
devote a separate paper to it, see [8]. 



2.3 Third model 



The third model is another one-parameter extension of the model treated in Subsection 12.11 
(different from that treated in the preceding Subsection I2.2p . Again its discrete-time equations 
of motion can be presented in three equivalent versions. 

The first is characterized by this prescription: the twice- updated N coordinates z n = z n (£+2) 
are the N roots of the following equation in the variable z, 



N 

E 

k=l 



z k - a + z k \ 



/v 



z - a + z k 



■) n 



Zk Zj 



1 



(2.12a) 



amounting again to the determination of the roots of a polynomial of degree N in the vari- 
able z. Here and below a+ and a_ are 2 arbitrary constants. 

Remark 2.6. As entailed by a comparison of these discrete-time second-order equations of 
motion with those of the first model, see (|2.1ap . this third model coincides, for a_ = 1, with the 
first model with a = a + . 
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A neater formulation of these equations of motion reads (after multiplication by a+, via (|A,10|) 
with Ck — r/j. = a + 5fc) as follows: 



N , 2 

z — a + Zj \ a_ + a + 



n 



z — a+Zj 



(2.12b) 



An equivalent, second formulation of this model is provided by the following system of TV 
polynomial equations for the twice-updated coordinates z n = z n (l + 2): 



N 

£ 

k=l 



Zk ~ a+Zk 



N 

n 



/ a + z k - Zj 



a_a. 



n = l,...,N; 



(2.13a) 



and a neater version of these equations of motion reads (again via (|A.10p . but now with Ck = z k, 
r]k = a+z k an d z replaced by a?\z n ) 



N 

n 



(a+ + a_)a 



N-l 

+ ' 



n = 1 N. 



\ Zj a^-z n 

3=1 \ J ^ / 

And a third, also equivalent, version of these equations of motion reads as follows: 

A' 



(2.13b) 



n 



-a_a+ -1 , n = l,...,N. 



(2.14) 



Remark 2.7. Remark 12.11 also holds for this model. 



The solution of the initial-value problem for this discrete-time dynamical system is provided 
by the following 

Proposition 2.3. The N coordinates z n (£) are the N eigenvalues of the N x N matrix 

U(£) = {a + fC+ + (a_)*C_, (2.15a) 

where the two constant {i.e., I -independent) N x N matrices C+ and C_ are defined in terms 
of the 2N initial data z n (0) and 5(0) = z n (l) by the formula 



C± = ±(o+ - a_)~ 1 [?7(l) - a T U(0)} 
with the two matrices £7(0) and £7(1) defined componentwise as follows: 



[£7(0)] nm = S nm z n (0), 

N 

[£7(l)] nm = a\~ N ' 



(2.15b) 



(2.15c) 



N 



^z fc (l)-a+ ^2 z fc (0) 

k=l k=l, k^n 



N 



* n 

j=l t j^m 



a+z m (0) - 2/(1) 



n, m = 1, . . . , N. 



(2.15d) 



z m {0)- Zj (0) 

Note that we are, for simplicity, assuming that the two coupling constants a± are different, 
^ a_ (see (j2l5b|) ). 
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It is plain from these formulas that, if the two "coupling constants" a± (are different and) 
are conveniently written as follows, 



with q± and r± real, then, if the two numbers q± are both nonnegative, q± > 0, the time 
evolution of the N x N matrix U{£) is bounded for all values of the discrete-time independent 
variable £ = 0,1, 2,..., hence its N eigenvalues z n {£) are all as well bounded (the motion is 
confined); if in particular the two numbers q± both vanish, q± = 0, and the two numbers r± 
are both rational numbers, r± = K±/L± with K + , L + and L_ coprime integers (and, 
for definiteness, L± > 0), then the discrete-time evolution of the matrix U(£) is periodic (with 
a period L independent of the initial data, being the minimum common multiple of L + and L_, 
L = mcm[L + ,L_]) hence the discrete-time dynamical system f|2. 12|) is isochronous; while if, of 
the two numbers q±, one vanishes and the other is positive, q+ = 0, q~ = q > respectively 
q— = 0, q + = q > 0, and r + respectively r_ are rational numbers, then the isochronous behavior 
(with period L + respectively L_) only emerges asymptotically, as £ — > oo, up to corrections of 
order exp(— q£). While clearly if q + and g_ are both positive entailing (see (|2.16p ) \a±\ < 0, 
then (see ()2.15ap ) the matrix U(£), hence as well all it eigenvalues z n (t), vanish asymptotically 
(as t — > oo): z n (oo) = 0, n = 1, . . . , N. 

Finally let us mention the transition from this discrete-time model to its continuous-time 
counterpart. The treatment is completely analogous to that detailed at the end of Subsection l2.lt 
except that now (|2,8bj) must be replaced by 

a± => 1 — iuj±e 



It is then easily seen that again, at order e° = 1, one gets from (|2.12b|) or (I2.13bp or (|2.14p a trivial 
identity, while at order e one gets the continuous-time goldfish equations of motion (|1.2aP : and 
the solution of this model, see Proposition 12. 3[ reproduces in this continuous-time limit the 
prescription ()1.2d[) . 

2.4 Fourth model 

The fourth model is also characterized by three equivalent versions of its discrete-time equations 
of motion. The first consists of the following prescription: the twice-updated N coordinates 
z n = z n(£ + 2) are the ./V roots of the following equation in the variable z, 



a± = exp(—q± + 2irir±) 



(2.16) 




with 



U! + = OtUJ 



0J- = 



(a — 




(2.17a) 



where the TV" quantities gk{z,z) are defined as follows: 





n = l,...,N. 



(2.17b) 
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Throughout this subsection, the following Subsection 13.41 and Appendix [Bl 

. VP , Pp 

a = a + - , b = - , 

1 — 7 1 — 7 

entailing 

a/3 = a/3-bq, rjp = (a - a)(l - 7), /3p = b(l--y), 



(2.17c) 



(2.17d) 



where a, (3, 7, rj and p are 5 arbitrary constants (but the 3 constants (3, rj, p only enter as f3p 
and ijp, hence any one of these three constants could be replaced by unity without significant 
loss of generality) ; in the following we use interchangeably these constants in order to simplify 
some formulas. 

An equivalent formulation of these discrete-time equations of motion reads as follows: 



N 

E 

k=l 



9k{z,z) 



(rjzk +j$)(z k - az n - b) 



(VZn + P) ' 



n = l,...,N, 



(2.18a) 



with 



9n{l,l) 



A-N 



7 



, Zn CLZn 



N 

") n 



n = l,...,N. 



(2.18b) 



It is a matter of trivial algebra to rewrite these equations of motion, (|2.18p . as follows: 

N 

1 



N 



E 

k=i 



I] {zj - azk ~ b) 



{rjz k + /3) (z k - az n - b) N , 

I! {Zj ~ z k) 

j=l, d+k 



7a 



JV-1 



r)z n + f3 ' 



n = l,...,N. (2.19a) 



And, as shown at the end of Appendix [Bl a neater version of this system of equations of motion 
then reads as follows: 



N 

n 



zj — a?z n — 6(1 + a) 



N 

n 

3=1 



rjZj + aj3 



7a 



N _ x riaz n + /3 + r]b 
r\z n + (3 



r}Zk+P / 

n = l,...,N. (2.19b) 
And a third, equivalent formulation of these equations of motion reads as follows: 

9n{l,z) = g n (z,z), n = l,...,N, (2.20) 
with g n (z, z) respectively g n (z,z) defined by (|2.18b|) respectively (|2.17b|) . 

Remark 2.8. Above and below we assume for simplicity that the parameters characterizing 
this model have generic values, for instance 7^0 and 7 7^ 1 (see (j2.17aj> and (|2,17cj) ) and a/0 
(see (|2l8b|) 1. 

Remark 2.9. This model - as the original goldfish model ()1.2ap . and as the models treated 
above, see Remarks 12.11 12.51 and !2.6l is invariant under a rescaling of the dependent variables, 
z n =>• cz n with c an arbitrary constant; but only provided the parameter (3 - hence as well the 
parameter b, see (|2.17cp - is also rescaled, j3 c/3, b =^ cb. 

The solution of the initial-value problem for this discrete-time dynamical system is provided 
by the following 
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Proposition 2.4. The N coordinates z n {€) are the N eigenvalues of the N x N matrix 

t 



U{£) = C/(0)P(0, 1 - 1) + [(B^' 1 + b ) P (K t-\% 



(2.21a) 



k=l 



where the N x N matrix P{1\^2) is defined as follows 



3=h 



(2.21b) 



and the two £-independent N x N matrices A and B are defined as follows 
A = iyV(0) - ^b, B = f3V(0) - b. 



(2.21c) 



Here and throughout we use the convention that {for arbitrary finite Xj) [ Xj = I if l\ > £2 



3=*1 



and Yl Xj = if k\ > hi- As for the two N x N matrices U(0) and V^O), they are defined in 

k=ki 

terms of the 2N initial data z n (0) and z n (0) = z n (l) as follows: 
U(0) = Z(0) =diag[z n (0)], 

V(0) = foZ(O) +/3]- 1 {M(0)Z(l)[Af(0)]- 1 -aZ(0)}, 

with the N x N matrices Z(£) and M(0) defined, componentwise, as follows: 

Z{£) = diag[z n (^)]; Z nm (£) = 5 nm z n (£), n,m = l,...,N, 
ffm(0) 



M nm (0) 



z m (l) - az n (0) - b 



n,m = 1, . . . , N, 



where the notation g m (0) is an abbreviation for g m (z,z), see (I2.17bj) . evaluated at z = z(0), 
z = 1(0) = z(l). Note that Lemma\AM (with f n = 1, g m = g m (0), £ m = z m (l), rj n = az n (0) + b) 
entails the following componentwise definition of the inverse matrix [M(0)]~ : 



{mo)]- 1 } nm = a 



l-N 



z n (l) - az m (0) - b 

9n{0) 



N 



n 

J=l, 3^n 



Zj(l) - az m (0) - b 



x < 



N 



n 



z n (l) - azj(0) - b 



Zj(l) - z n (l) 
> , n, m = 1, . . . , N. 



z m (0)- Zj (0) 

hence an explicit expression of the N x TV matrix V(0) reads, componentwise, as follows: 



V ( 0) - az "(°) 5 + 

v nrn\ {J J — /n \ , n u nm 



l-N 



N 

V z n (0)+P~'"' L ' r/z n (0)+/3^ 



( N 

II [z k (l)-a Zj (0)-b] 

j=l, ii=n 



N 

**o){ n 

3=1, 3+k 



Zj(l) - az m (0) - b 
Zj (l)-z k (l) 



x < 



n [z m (o) - Zj (o)] 

3=1, j¥=™ 
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Remark 2.10. It is relevant to this expression, (|2,2ip . of the N x N matrix U(£) - whose N 
eigenvalues provide the N coordinates z n (£) - that (|2.21b|) and (|2.21c|) entail 

P(£i,£ 2 ) = Qdia,g[ Pn (£ 1 ,£ 2 )]Q- 1 , A -- 

£2 

Pn(£i,£ 2 ) = JJ {a n n j + a), a n = r,v n 
Hi 

with v n the N (^-independent) eigenvalues of the N x N matrix V(0) and Q the corresponding 
(^-independent) diagonalizing matrix, 

V(0) = Qdmg(v n )Q-\ 

And let us mention that, also for this fourth model, a transition from the discrete-time 
independent variable £ to the continuous-time variable t can be performed (see the end of 
Subsection I3.4p . And, as in the case of the second model, also in this case the continuous-time 
goldfish model thereby obtained turned out to be, to the best of our knowledge, new. Hence it 
seemed appropriate to devote to this model a separate paper [9]. 



= Qdiag[a n ]Q , B = Q dia,g[b n ]Q , 
V 

--6, b n = (3v n -b, n = l,...,N, 



3 Proofs 

In this section we prove the findings reported in the preceding Section [2j 

The basic strategy to obtain all these results goes as follows. The starting point is a solvable 
system of two matrix first-order discrete-time ODEs, say 

U = F 1 (U,V), V = F 2 (U,V), (3.1) 

where £ = 0, 1,2, ... is the discrete-time independent variable, the two dependent variables 
U = U(£), V = V{£) are N x N matrices and of course superimposed tildes denote the updating 
of the discrete-time, U = U(£ +1), V = V{£ + 1). The solvable character of this matrix system 
entails the possibility to obtain explicitly the solution of its initial-value problem. Four cases 
when this is possible - corresponding to 4 simple assignments of the functions Fi(U,V) and 
F2QJ, V) - are treated in the following 4 subsections. Note that the two functions Fi(U,V), 
F 2 (U, V) are assumed to depend on no other matrix besides U and V (and the unit matrix I); 
they may of course feature some scalar constants, and the order in which the two, generally 
noncommuting, matrices U and V appear in their definition is of course relevant: see below. 

One assumes moreover that the N x N matrix U = U(£) is diagonalizable and denotes as 
R = R{£) the diagonalizing N x N matrix: 

U^RZR' 1 , U{£) = R{£)Z{£)[R{£)]- 1 , (3.2a) 
Z = diag[z n ], Z{£) = dmg[z n (£)l (3.2b) 

where the notation z n {£) for the N eigenvalues of the NxN matrix U = U(£) shall be justified by 
the identification, see below, of these quantities with the dependent variables of the discrete-time 
dynamical systems introduced above. 

Remark 3.1. These formulas entail that the matrix R(£) is defined up to right-multiplication 
by an arbitrary diagonal matrix D(£), R{£) =>• R(£)D(£). 

Next we introduce the two matrices M(£) and Y(£) defined as follows: 



M = R- X R, M(£) = [R(£)]- l R(£ + 1), 



(3.3a) 
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V = RYR-\ V(£) = R(£)Y(£)[R(£ + l)]-\ (3.3b) 
so that 

V = RYM~ 1 R- 1 , V{£) = R{£)Y{£)[M{£)]- l [R{£)Y l . (3.3c) 



Remark 3.2. The element of freedom in the definition of the matrix R{£), see Remark 13.11 
entails that the matrix M{£) is defined up to the "gauge transformation" resulting by inserting 
in its definition (pOa]) the N x N matrix [R(e)D(£)]- 1 = [£>(f)] _1 [R^)]' 1 in place of the matrix 
[i?(^)]~ 1 (and of course R(£ + 1)D(£ + 1) in place of R(£ + 1)): hence, as a consequence of the 
arbitrary nature of the diagonal matrix D(£), out of the N 2 elements of the N x N matrix 
M = M(£) only N 2 — N are significant. Likewise for the matrix Y = Y(£). 

One then, by inserting (|3.2a|) and (|3.3c|) in (|3.1|) . obtains the following system of two first- 
order discrete-time N x N matrix evolution equations: 

MZ = Fi{Z,YM~ l )M, MY = F 2 (Z,YM~ 1 )MM; (3.4) 

and from these two matrix equations, by making a convenient ansatz for the two matrices M = 
M{£) and Y = Y{£) in terms of the 2N quantities z n = z n {£) and z n = z n {£ + 1) - an ansatz 
which must of course be consistent with these two matrix evolution equations - one obtains 
a system of N second-order discrete-time evolution equations for the N coordinates z n = z n (£). 
This last step is of course only possible for special assignments, in the discrete-time matrix 
evolution equations (13. lj) . of the two matrix functions F\(U,V) and ^(U, V), see below. 

The discrete-time dynamical system thereby obtained is then solvable, since the quantities 
(£) are the N eigenvalues of the N x N matrix U = U{£) which, as solution of the, 
assumedly solvable, matrix evolution system (|3.ip . can be explicitly evaluated. How this works 
out is shown in detail in the following subsections: in more detail in Subsection 13.11 where the 
simplest case is treated. 

Let us also mention, once and for all, that in the following we will conveniently assume that 
the matrix U is initially diagonal: 

U(0) =Z(0) = diag[z n (0)], (3.5a) 
implying (up to the ambiguity mentioned above, see Remark 13. ip 

R{0) = I. (3.5b) 
Here and throughout / is the N x N unit matrix, i.e., componentwise, I nm = 5 nm . 

3.1 Solution of the first model 

The point of departure to obtain the findings reported in Subsection 12. II is the following discrete- 
time first-order, linear, matrix system (see (13.ip ): 

U = aU + V, V = V, (3.6a) 

where a is an arbitrary scalar constant. Note that the second of these two ODEs entails that in 
this case V is a constant (i.e., ^-independent) N x N matrix, V{£) = V(0). It is plain that the 
solution of the corresponding initial-value problem for the N x N matrix U reads 

n l - 1 

U(£) = U(0)a e + V(0)- -. (3.6b) 
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Let us now proceed as indicated in the first part of Section O It is then easily seen 
(via (|3.2p and (|3.3|) ) that the first of the two discrete-time matrix evolution equations (|3.6a| 
yields (see (|3.4p ) the matrix equation 

MZ - aZM = Y, (3.7a) 

namely, componentwise, 

M nm = „ Y ™ , n,m = l,...,N. (3.7b) 

Likewise, the second of the two discrete-time matrix evolution equations (|3.6ap yields the 
matrix relation 

YM = MY. 

Via (|3,7b|) this matrix equation implies the following iV 2 relations: 



JV 



I Y nk Y kr 

k=l 



az k ) -{z k -az n ) 1=0, n, m = 1, . . . , N. (3. 



This derivation shows that this system of N 2 discrete-time equations of motion is equivalent 
to the solvable equation of motion (|3.6ap for the NxN matrix U; hence it is just as solvable. Note 
that the dependent variables are now the N coordinates z n and the N 2 matrix elements Y nrn 
(of which only N(N — 1) are significant, see Remark 13.2} so the number of equations and the 
number of dependent variables tally). To obtain a model that qualifies as discrete-time analog 
of the continuous-time goldfish model (]2.9aP we need to distill from this system a set of only N 
equations of motion involving only the N coordinates z n . The standard trick to do so (see, for 
instance, Section 4.2.2 entitled "Goldfishing" of [6]) is to identify - if possible - an ansatz which 
expresses the A^ 2 components of the matrix Y in terms of the 2iV quantities z n , z n , yielding 
iV equations of motion involving only the N coordinates z n , z n and z n - to be interpreted 
as equations of motion of the discrete-time goldfish - and implying that the N 2 equations of 
motion (|3.8p are all satisfied, thanks to these very equations of motion. 

An educated guess for such an ansatz reads as follows: 

Y nm = gm, n,m = 1, . . . , N. (3.9) 

Note that we reserve at this stage the option to assign the quantities g m . 
Via this ansatz the equations (|3.8|) become 



jv , 

yl^Jh _Jh ) =0 , n,m = l,...,N, (3.10a) 



k=l 



z m - az k z k ~ az n 



hence they amount to the following 2 systems, each involving only N equations: 

N , . 

V L 9k = 1, n = l,...,N, (3.10b) 
^ \ z k - az n I 



k=l 

N 



k=l 



E M^T H 1 ' n = l,...,N. (3.10c) 

x z n - az k J 



The unit in the right-hand sides could of course be replaced by an arbitrary constant c - of 
course the same constant in (|3.10bp and (|3.10cp - but this would merely entail an irrelevant 
rescaling of g k by c; see below. 
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These are now two sets of N equations, each featuring linearly the iV quantities gk, that 
we like to eliminate in order to obtain a set of N "equations of motion" determining the twice 
updated coordinates z n = z n (£+2) in terms of the 2N coordinates z n = z n {£) and z n = z n (£+l). 
There are three alternative strategies to achieve this goal. One can solve the first linear system 
thereby obtaining g k as a function of z and z, and then insert this expression g k = gk (z, z) in 
the second system; alternatively, one can solve the second linear system, thereby obtaining gk 
as a function of z and 2, and then insert this expression gk = gk{z_, z) in the first system; or 
one can equate the two expressions of gk obtained solving the first, respectively the second, 
system, i.e. write gk(z,z) = gk(z,z). Clearly the three sets of equations of motion obtained 
in this manner are equivalent, i.e. they characterize the same discrete-time evolution of the N 
coordinates z n = z n (£); but they may seem quite different (indeed, see (|2.1|) . (|2.2|) and (|2.3|) ). 
Note that we introduced a superimposed decoration on the functions §k(z, z) respectively gk (z, z) 
to emphasize that the functional dependence on their arguments is generally different, as implied 
by their definitions as solutions of (|3. 10b[) respectively of (|3.10cp . 

Let us first of all see what the first approach yields. From (I3.10bj) one obtains (via Lemma [A. II 
reported in Appendix[A] with = z k , f] n = az n , c = 1) the following expression of g k = g k (z,z): 

9k(z,l) = (zk-az k ) I ~ k = l,...,N. (3.11) 

The insertion of this expression of gk in (j3. 10c[) yields the equations of motions (|2,lap . 

Likewise, the second approach yields, from (|3. 10c|) (again via Lemma IA.11 but now with 
£ fc = az k , r} n = z n , c = -1) the following expression of g k = gk{z,z): 



g k (z,z)=a 1 ~ N (zk-azk) ^— =^ , k = l,...,N. (3.12) 




The second version, (|2.2ap . of the discrete-time equations of motion follows by inserting this 
expression of g k in (|3.10b|) . 

And the third approach yields, by equating (|3.1ip to (|3.12p . the third version, (|2.3p . of the 
discrete-time equations of motion. 

We have seen that the solutions z n (£) of these discrete-time equations of motion are provided 
by the eigenvalues of the N x N matrix U {£) , see (j3.6b[) . To prove Proposition 12.11 we must 
now obtain from (|3.6bj) (also taking advantage of the ansatz (|3.9p ) the expression (|2.4p of this 
matrix in terms of the initial data z n (0), z n (0) = z n (l) of the discrete-time dynamical system. 

This requires that we express the two matrices U(0) and V(0) appearing in the right-hand 
side of (I3.6bj) in terms of the initial data z n (0), z n (0) = z n {l). 

The expression of U(0) is an immediate consequence of f)3.5a[) : 

[U(0)]nm = S nm Zn(0). (3.13) 

To obtain V(0) we note first of all that (|3.3cp and (|3.5bj) imply 
V(0) = Y(0)[M(0)]-\ 
while (I3.7bp with the ansatz (13. 9h implies (at £ = 0) 

M nm (0)= 9m{0) n,m = l,...,N, 

z m (l) - az n (0) 

where of course g m (0) stands for g m (z,z), see ([3. lip , evaluated at z = z(0), z = z(0) = z(l). 
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We then evaluate the matrix [M(O)]" 1 via Lemma IA.5I (with £ m = z m (l), rj n = a z n (0), 
f n = l and g m = g m (0)) and, using again the ansatz (|3.9p (at £ = 0), we obtain the following 
expression of the N x N matrix V(0): 



[V(0)] nm = v m (0)u nm , n,m = l,...,N, 
with v m (0) defined as in Subsection l2,ll (see (|2.4bp and the sentence following this formula) and 



N 



II,, 



n [z k (l) - a Zj (0)] 

3=1, i+m 



N 

£ 

1 n fed) - 

3=1, j¥=k 



N 



n, m = 1, . . . , N. 



But the identity (|A.8p (with r)k = £fe(l), Q = az j(fi)) entails u nm = 1, hence 
[V(0)] nm = v m (0), n,m = l,...,N. 



(3.14) 



The insertion of these expressions of U(0) and V(0), (|37L3j) and ([331]), in (|3T6b|) yields (p^]) . 
thereby completing the proof of Proposition 12.11 

Let us now provide a terse treatment of the transition from the discrete-time equations of 
motion (|2.1b|) . which we conveniently re- write here as follows, 



JV 

n 



1 + a, 



n=l,...,N, 



(3.15) 



a Zn 



to the continuous-time case, see (I2.9ap . It is then appropriate to treat separately the factor with 
j = n in the product appearing in the left-hand side of (|3.15p . and all the other factors with 
j ^ n. The basic equations are (|2,8p . entailing 



a 2 z 



n = 2e{z n + iujz n ) + e 2 (2z n +uj 2 z n ) + 0(e 3 ), 



z n - az n = e(z n + iuJZn) + — (3'z n + 2iuz n ) + 0(e 3 ), 
i 2 zj = z n — zj + 2e{z n + iuiZj) + 0(e 2 ) , j ^ 



z n - azj = z n - zj + e(2z n - zj + iuzj) + 0(e 2 ) , j ^ n. 
Hence, after a little algebra, 

, - 2iujz n + uj 2 z ri 



a 2 z n 



Zn Q/Z r , 



a Zj 



Zri, CLZ-j 



--2 + e 
l + e 



Z n + iL0Z n 



+ 0(e 2 ), 



Zn Zj 



implying 

N 



n 

3=1 









/ Z n 


-a 2 Zj \ 






- azj ) 


-( 



2 + e- 



+ 0(e*), j^n, 



— 2iuz n + uj 2 z n 

Z n + iiOZ n 
N 



(3.16a) 
(3.16b) 



3=1, 3 5^1 



+ lUZj 
Zn Zj 



+ 0{e 2 ). 



(3.16c) 
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While of course 

1 + a = 2 — iue, 



(3.16d) 



see ()2.8b|) . It is then clear that the insertion of these two formulas, (|3. 16cj) and (|3.16d|) . in (|3.15p 
yields, at order e° = 1, the trivial identity 2 = 2, and at order e the equations of motion of the 
continuous-time goldfish model (|2.9ap . 

In an analogous manner one reobtains (I2.9ap from (l2.2bp or from (12.3p . 

Let us also show that (|2.5p . which we rewrite here conveniently as follows, 



N 

n 

k=l 



z - a £ z k (0) 



a a 



1 



(3.17) 



yields, in the continuous-time limit, (|2.9b|) . Indeed the relation a = 1 — ieuj (see (|2.8bj) ) entails 
cr 1 = 1 + ieuj + 0(e 2 ), 



exp(-iut) I 1 — e 



T 



+ o(. 2 ) 



(via the first of the three relations (|2.8ap ). and 

z k (l) = z k (0)+ez(0) + O{e 2 ) 
(via the second of the three relations (|2.8cp . with £ = 0). Via these three relations (|3,17p becomes 
Y z - exp(-i^)(l - eu 2 t/2){z k (0) + e[z k {0) + iuz k {0)}} + 0(e 2 )' 



n 

k=l 



z - exp(-iwt) (1 - eu 2 t/2)z k (0) + 0(e 2 ) 
exp(-iwt)(l - euj 2 t/2){l + ieuj) - 1 + 0(e 2 ) 



exp(-iujt) (1 - euH/2) -l + O (e 2 ) 

i.e. (dividing each numerator by the corresponding denominator) 

N 



n 

k=l 



1-e fe(0)+^(0)] +o(g2) 
z — exp(— iu)t)z k (0) 



1 + e 



iuj 



ex.p(—iut) 



J + Oie 2 ). 



Clearly to order e° = 1 this yields the trivial identity 1 = 1, and to order e just the formula (I2.9bp . 

Let us end this subsection by pointing out that there is another ansatz that allows to trans- 
form the system of N 2 equations (I3.8P into two separate systems of N equations, but only in 
the special case a = 1. This alternative ansatz reads (instead of (|3.9p ) 



Y 



fn 



n,m = 1, . . . , N, 



entailing (but only provided a = 1) the replacement of the system of TV 2 equations (|3,8p with 
the following two systems of N equations: 



N 



E 

k=l 

N 

E 

k=l 



fk 



(Zk ~ Z n ) 2 
fk 



z k y 



1, n = l,...,N, 



1, n = l,...,N. 



But, as indicated at the end of Section [TJ we postpone the treatment of the corresponding class 
of discrete-time dynamical systems to a separate paper. 
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3.2 Solution of the second model 



The proof of the findings reported in Subsection 12.21 is analogous to that provided above, see 
Subsection 13.14 so our treatment in this subsection is quite terse, being limited to indicate the 
changes with respect to that reported in the preceding Subsection [3TTJ Now the system of matrix 
evolution equations (I3.6a|) is generalized to read 



U = U(al + bV) + V, V = V; 



(3.18) 



hence its solution is given by (|2,llap . Clearly this evolution equation, (|3.18p . respectively its 
solution, (12. Haft , reduce to (|3.6aj) respectively to (I3.6bj) when b vanishes. 
The rest of the treatment is analogous. (|3.7ap is now generalized to read 

MZ - aZM = (J + bZ)Y, 

hence it yields, in place of (|3.7bj) . 

1 + bZr, 



Mr, 



Y 



n, m 



1. 



,N. 



, QZ n 

In place of (|3. lOaj) (again via the ansatz (|3.9p ) one now has 

N 



E 

k=l 



g k (l + bz k ) g k (l + bz n ) 



az k ) 



z k - az % 



0. 



n, m 



1,...,N, 



(3.19a) 



hence in place of ()3.10bp and (|3. lQcj) one gets the two sets of N equations 

N 



E 

k=l 

N 

E 

k=l 



9k 



1 



g k (l + bz k 



az k 



l + bz n ' 
= h 



n 



n 



N. 



(3.19b) 
(3.19c) 



By solving the first set one obtains (via Lemma [A. 21 with £ k = z k , rj n = az n , c n = l/(l + bz n ), 
and then the identity ()A.7j) with z = —a/b, rj k = az k , (j = Zj) the following expression of 
9k = 9kU,z): 



9k{z,z) 



Zk ~ azk 
1 + bz k 



N 

n 



Zk - azj 

Z k Zj 



1 + bzj/a 
1 + bzj 



(3.20a) 



By solving instead the second set one obtains (via Lemma lA.H with £ k = —az k , r] n 
c = 1 and g k replaced by g k (l + bz k )) the following expression of g k = g k (z,z): 



9k (I, I) 



1— JV / z k ~ az k 




(3.20b) 



The three versions, (|2.10p . of the equations of motion reported in Subsection 12.21 then follow 
by inserting (I3.20al) in (|3.19cj> . by inserting (I3.20bl) in (I3.19bl) . and by equating (I3.20al) to (|3.20bj) . 

Next, let us prove Proposition 12.21 One proceeds again in close analogy to the treatment 
of the preceding Subsection 13.11 hence we only mention where the treatment here differs from 
that provided there. It is easily seen that the expression of ?7(0) is the same as that given 
there, see (|3.13p . while the expression of V(0) (because one must now use Lemma [A.5I with 
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f n = 1 + bz n (0) rather than f n = 1) is now given by ()2. 1 lc|) . The insertion of these expressions 
of £/(0) and ^(0) in (|3.18|) reproduce (|2.11|) . thereby proving Proposition 12.21 

Let us end this subsection by outlining what happens in the continuous-time limit which 
obtains by setting 

a = 1 + ebrj, V(0) = eB 

with e infinitesimal, and correspondingly replacing the discrete-time matrix evolution equa- 
tion ([37T8]) with the matrix ODE 

U = bU(r]I + B) + B, (3.21a) 

the solution of which reads 

U(t) = U(0) exp[b(T]I + B)t] + B[b(rjl + S)]~ 1 {exp[6(? ? / + B)t] - I}. (3.21b) 

As already mentioned in Subsection 12.21 the (continuous-time) goldfish-type model obtainable 
by focussing appropriately on the evolution of the N eigenvalues z n (t) of this N x N matrix U (t) 
(evolving according to (|3.21aj) ) was, to the best of our knowledge, new, when the solvable matrix 
evolution equation (|3.21a|) was identified as continuous-time limit of (|3.18|) ; its treatment is 
provided in [8]. 



3.3 Solution of the third model 

The starting point is the following linear system of two discrete-time matrix evolution equations: 

U = a + U + PV, V = a-V, (3.22a) 

where the 3 constants a±, f3 are a priori arbitrary (/5 ^ 0). It is easily seen that the solution of 
the initial- value problem for U reads as follows (with an analogous formula for V): 

U(£) = a l + C + + a e _C_, (3.22b) 

and the two constant matrices C± given by (|2.15bj) . 

Remark 3.3. This solution U (I) of the initial- value problem for the discrete-time N x N 
matrix evolution equation (|3.22ap depends only on the 2 constants a±: see ()3.22bp and (j2. 15b[) . 
Indeed the system of two first-order discrete-time evolution equations (13.22aj) is easily seen to 
correspond to the single second-order evolution equation 

U - (o + + a_)C7 + a + a_U = 0, 

from which the constant (3 has disappeared (but note that this second-order matrix ODE obtains 
only if f3 / 0; indeed if /3 = 0, U satisfies a first-order evolution equation, see the first of the 
two equations (|3.22a[) . 

Remark 3.4. For a+ = a, j3 = 1, a_ = 1, the system (|3.22ap coincides with (|3.6ap hence this 
model reduces to the first model, confirming Remark 12.61 

We then proceed as in the first part of Section [3l It is then easily seen that the two matrix 
evolution equations (|3.22ap become 



MZ = a + ZM + PY, 
MY = a-YM. 
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Via (|3.2bp the first of these two equations reads, componentwise, 

M nm = . PYnm , n,m = l,...,N, (3.23) 

Zm 0'+Z n 

and using this formula it is easily seen that the second can be written, componentwise, as follows: 



k=l 



Y nk Y k 

i. 



1 



z m - a + z k z k ~ a + z r, 



0, n,m = l,...,N. (3.24) 



At this point we use again the ansatz (|3.9[) for the matrix Y nm , the consistency of which 
is vindicated by the subsequent developments. Here the N quantities g m are again a pri- 
ori arbitrary; they shall be determined as functions of the 2N un-updated and once-updated 
coordinates z n = z n {€) and z n = z n (£ + 1), or alternatively of the once and twice updated coor- 
dinates z n = z n (£ + 1) and z n = z n (£ + 2), see below. It is indeed immediately seen that via the 
ansatz (|3.9|) the system of N 2 equations (I3.24|) becomes 

N / \ 

Y^(^a-9k _^_9k =Qj n>m = i j ...,jv ; (3.25a) 

k=1 \z m -a + z k z k -a + z n J 

hence it can be replaced by the following two separated systems of only N equations: 

E =1 ' n = l,...,N, (3.25b) 

fr[ \ z k - a + z n J 

N ( \ 1 

E U - ) = — n = l,...,N. (3.25c) 

^[\z n -a + z k / a- 

The first, ()3.25b|) . of these two systems defines uniquely the quantities g k = g k (z_,z), 
yielding again, via (|A.13jl . the expression (13. lip (with a replaced by a+): 

N f~ \ 

9k(z,z) = (z k -a + z k ) H Z \ _t 4 , k = l,...,N. (3.26) 

3=l,^fc V ^ Zj J 

Insertion of this expression in the second, (|3.25c|) . of the two systems written just above then 
yields the evolution equation (I2.12ap . The identification of the third discrete-time dynamical 
system of goldfish type, see ()2.12ap . is thereby accomplished. 

The second version, (|2.13ap . of this model obtains by solving for g k = g k (z_,z) the se- 
cond, (|3.25cp . of the two systems written above (via Lemma I A. 11 with £ k = a+z k , r] n = z n 
and c = — l/a_), thereby obtaining 

_ (LzSta) n ( a ^~\ ), * = !,...,»; (3.27) 

and by then inserting this expression of g k in (|3.25bj) . 

And the third version, (|2.14p . of this model obtains by equating (|3.26p to (|3.27p . 

As for the proof of Proposition 12.31 it follows immediately from the above treatment, see 
in particular (|2.15p : there only remains to justify the identification of the two matrices U(0) 
and U(l), see fl2TT5c]) and (|2lM|> . 

The first of the two formulas, (|2.15cj) . is just (|3.5a|) . 
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To prove the second, (j2.15d|) . we note that (|3.5ap and (|3.3ap entail (for t = 0) 

12(1) = Af(0), (3.28) 

while ([32]) (with £ = 0) reads 

17(1) = R(1)Z(1)[R(1)]- 1 . (3.29a) 

Hence (via (13381) ) 

17(1) =M(0)Z(1)[M(0)] -1 . (3.29b) 

We now note that (|3.23[) and the ansatz (|3.9p imply that the N x N matrix M(0) is defined 
componentwise as follows: 



M nm (0) 



9m(0) 



z m (l) - a + z n (0)' 



n,m = 1, . . . , N, 



(3.30a) 



with (see (137261) ') 



N 



g m (0) = [Zm(X) - a+z m (0)] 



so that 



N 



z m (l) - a + z n (0) . 



z m (l) - a+Zj(0) 

Z m (l) - Zj(l) 



z m (l) - a + Zj(0) 

Z m (l) - Zj(l) 



m = l,...,N, 



(3.30b) 



n,m = l,...,N. (3.30c) 



Next, we note that the expression (|3.30ap of the matrix M(0) entails - via Lemma [A. 51 (with 
fn = 1, 9m = 5m(0), Cm = 2m (1), ??n = a+^n(°)) and th e expression of g m (0) given above - that 
its inverse, appearing in the right-hand-side of (|3.29bp . is explicitly given, componentwise, as 
follows: 



{l M (0T 1 } nm = a 
n, m = 1, . . . , N 



l-N 
+ 



0,+ Zm (0) - z m (l) 
(0) - z n (l) 



N 



n 



a + z m (0) - Zj(l) 



z m (0)- Zj (0) 



(3.30d) 



The insertion of this formula and (|3.30c|) in (j3.29bp entails that the matrix U(l) reads compo- 
nentwise as follows: 



N 



U nm (l) = al~ N 



i=l, jj^m 
( 



Q>+Zm 

(O)-^(l) 
z m (0) - Zj (0) 



N 

E 

k=i 



z k (l) 



J! [z k (l)-a + Zj(0)} 

n N(i)-^-(i)] 



n, m = 1, . . . , N. 



And it is then immediately seen that this formula yields, via the identity (|A. 12j) (with % = Zfc(l) 
k = 1,. . . ,N; Q = a+Zj(0), j = 1, . . . ,n - l,n + 1, . . . ,N), the formula (j2.15d|) . 
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3.4 Solution of the fourth model 

The treatment in this subsection is rather terse, since it is analogous to that of the preceding 
subsections; and the notation is of course analogous. But now the starting point is the following 
nonlinear system of two discrete-time matrix evolution equations: 



U = aU + PV + vUV, V = p + 7 V, 



(3.31) 



featuring the 5 arbitrary constants a, (3, 77, p, 7 (which, as noted in Subsection 12.41 can De 
reduced to 4 by taking advantage of the freedom to rescale V) . 

It is a standard task to see that the solution of the initial-value problem for this matrix 
system reads as follows: 



V{£) = 7 V(0) + p 



7 



7-1 



with U{£) given by (22U>. 

We now proceed again as in the first part of Sectional Via (|3.3aj) and (j3.3c|) we get from (13. 31 1) 
the two matrix equations 



MZ = aZM + (J3 + r)Z)Y, 
MY = pMM + 7YM. 

The first, (|3.32ap . of these two matrix equations entails, componentwise, 



(3.32a) 
(3.32b) 



Y 



az r . 



r]z n + @ 



n, m 



hence the second, (|3.32b|) . of these two matrix equations yields (when written componentwise) 
the following N 2 equations: 



N 

E 

k=l 



M n kM km 



az k 



7]Z k +{} 



p-r 



Zk - az n 



0. 



n, m 



1 



)••• r'j 



N, 



(3.33a) 



which, as can be easily verified, can be conveniently rewritten as follows: 

N 



E 

k=l L 



M nk M km 



az k -b z k - az Tl 
-7- 



rjz k + /3 



r]Zn + P 



0. 



n,m = 1, . . . , N, 



(3.33b) 



with the two constants a and b defined by fl2.17c|) . 

Next, we make the following ansatz for the matrix M n 



Mr, 



9m 



n,m = 1, 



N. 



az r , 



(3.34) 



Here the N quantities g m are a priori arbitrary; they shall be determined as functions of the 
2iV un-updated and once-updated coordinates z n = z n (£) and z n = z n (£ + l), or of the once and 
twice updated coordinates z n = z n (£+l) and z n = z n (£ + 2), see below. It is indeed immediately 
seen that via this ansatz ()3.34p the system of N 2 equations (|3,33b|) can be rewritten as follows: 



N 

E 

k=l 



9kij]Zn + P) 



9k 



{rjz k +/3)(z k - az n - b) {z m - a z k -b) 



n,m = 1, . . . , N; 



(3.35a) 
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hence it can be replaced by the following two separated systems of only N equations: 

iV 



9k 



£ 

k=l 

N 

T 

k=1 \z n ~ az k -bj 7 



(r)z k + P)(z k - az n - b) 
9k 



n = l,...,N, (3.35b) 



{riz n + (3) ' 

n = l,...,N. (3.35c) 



The first of these two systems defines uniquely the N quantities g k = (jk(z,z), yielding their 
expression (|2.17b|) (see Appendix [B] for a proof). The second equation then yields the evolution 
equation (I2.17a|) . 

The alternative possibility is to determine (via Lemma [A. 11 with = a z k , rj n = z n — b, c = 
— 1/7) the quantities g k = g k (5, z) as solutions of the second system, yielding the formula (|2.18bj) ; 
and to then insert these expressions of g k = g k (z_, zj in the first system of equations. Clearly in 
this manner one arrives at the equations of motion (|2.19ap . 

And a third possibility is of course to equate (I2.17b[) to (12.18bl) . see ([230j) . 

The identification of the three variants of the equations of motion of the fourth discrete-time 
dynamical system of goldfish type, see Subsection 12.41 is thereby accomplished. 

The proof of Proposition 12.41 follows immediately from the above treatment; and we trust 
that the identification in terms of the 2N initial data z n (0) and z n {\) of the two matrices U(0) 
and V(0), see (|2.4j) and (|2,4p . is sufficiently obvious (also in the light of the analogous treatment 
in the preceding subsections of this section) not to require an explicit justification here. 

We end this subsection with a terse mention of the continuous-time model that obtains from 
that treated above (in this subsection) via the limiting transition from discrete to continuous 
time. The point of departure for the treatment of the continuous-time dynamical system of 
goldfish type obtained in this manner is the following continuous-time system of two first-order 
matrix evolution equations 

U = ai U + a 2 V + a 3 UV, V = a 4 + a 5 V, 

that obtains from (|3.3ip via the assignments t =^ el, U(£) U(t), V(£) =^> V(t), a = 1 + ea±, 
j3 = eo2, f] = £03, p = ea^, 7 = 1 + ea^, with e infinitesimal. The resulting model of goldfish 
type was, to the best of our knowledge, new] a detailed treatment of it is provided in [9]. 



4 Outlook 

In this paper we have introduced and tersely analyzed 4 different discrete-time dynamical 
systems of goldfish type. The possibility to identify other discrete-time evolution equations 
amenable to exact treatment by variations of the methodology used in this paper is open: let 
us outline here an avenue to such generalizations. 

Consider the system of two N x N matrix discrete-time first-order evolution equations 

Si S2 

[FiAU)F 2 AU)] = £ [*i,,{U)V$2s{U)], (4-la) 

S = l 8=1 

]T [F 3)S (U)F A)S (U)]V = $ 3 (U)V$ 4 (U), (4.1b) 

s=l 

where the two TV x N matrices U = U(£) and V(£) are the dependent variables, I = 0, 1, 2, . . . 
is the independent discrete-time variable, Si, S2, S3 are 3 arbitrary positive integers, iq ;S (ii), 
F2, s (u), F 3jS (u), F 4 , s (u) and $ M (u), $ 2 , fl (u), *s(«), *4(«) are 2(S 1 + S 2 + S 3 + 1) a priori 
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arbitrary (scalar) functions of their (scalar) argument u (of course becoming N x N matrices 
when the scalar u is replaced by an N x N matrix). Then introduce the eigenvalues z n (£) of the 
matrix U(£), as well as the matrices Z = Z(£), R = R(£), M = M{£) and Y = Y(£), as above 
(see (|3.2p and (|3.3p ). It is then plain that the matrix evolution equation (|4.1aj) becomes 



Si S 2 

]T [F ltB (Z)MF 2tS (Z)] = [^iAZ)Y^2,s{Z)] 

s=l 8=1 
entailing componentwise 



M, 



E [*i, s (%)$2, s (i 

s=l 



nm — g 

E [Fi, s {z n )F 2>s (z m )] 



Y 

± n. 



n,m= 1, . . . , N. 



8=1 



Likewise (|4.1bj) becomes 



S 3 



^ [F 3>s (Z)MF 4iS (Z))Y = <S> 3 (Z)Y$ 4 (Z)M, 



8=1 



entailing componentwise (via (|4.2bp ) 



(4.2a) 



(4.2b) 



N 



k=l 



S-2 



J2 F 3,s^n)F 4 , s (z k ) 



8 = 1 



IV 



$3(0 X/ 



fc=l 



YnkYkm'&iiZk) 



E [*l,a(^)*2,a(5fc)] 

cr=l 

Si 

E [i ? l,a(^)i ? 2 !f r(5fc)] 

<T=1 

E [$l,»(2k)*2,«(5m)] 



E [Fi >8 (z fc )F 3 , a (l m )] 



n, m = 1, . . . , N. 



And via the ansatz Y nm = g m (see (|3.9p ) this system of N 2 equations can clearly be replaced 
by the following two systems of N linear algebraic equations for the N quantities g k : 



N 

E 

k=i 



N 

E 

fc=i 



9k 



S 2 



^2F 3tS (z n )F iiS (z k ) 



.8=1 



S 2 

E [$lA Z n)$2,a(zk) 
o=l 

E [Fl,a(Zn)F2,a(Zk) 
o=l 



<!> 3 (z n ), n = 1,..., AT, (4.3a) 



S> 



E [*l, fl («fe)^2,s(«n)] 



8 = 1 



E [^1,8(4)^2,8^)] 



n = 1, 



iV. 



(4.3b) 



One can then solve the first, (|4.3ap . respectively the second, ()4.3bp . of these two linear systems 
for the N quantities g k , getting the expressions g k = gk(z,z_) respectively g k = gk{z_,£). One 
arrives thereby at three equivalent systems of N discrete-time second-order evolution equations 
for the N coordinates z n {£): (i) by inserting the expression gk(z,z) in (|4.3bp ; (ii) by inserting the 
expression gk(z,£) in (|4.3a|) ; (in) by setting g n (z,z) = g n (z,z^. While of course the evolution 
in the discrete time £ entailed by these equations of motions corresponds to the evolution of 
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the eigenvalues of the matrix U(£) solution of the matrix evolution equation (|4.1ap (with the 
ansatz (|3.9|) properly taken into account). Hence if that matrix evolution equation, (|4.1a|) . is 
solvable, the discrete-time dynamical system described by these three equivalent sets of second- 
order evolution equations is as well solvable. 

One has thereby identified a discrete-time solvable dynamical system. A remaining open 
question is the extent to which its equations of motion can be exhibited in reasonably neat 
form: this depends on the extent that the two quantities gk = gk {z, 1) respectively gk = gk (l> l) 
defined as solutions of the two linear systems (j4.3a|) respectively (|4.3b|) can be expressed more 
explicitly than via the standard Cramer formula (ratio of two determinants). 

Clearly the 4 models treated in this paper belong to this class (14. ip (the last one, however, 
only if p = 0): see (|3.6aj) . (|3.18p . (|3.22aj) and (|3.3ip . The interest of additional models treatable 
via this approach depends on the neatness of the corresponding equations of motion, which can 
only be investigated on a case-by-case basis. 



A Appendix 

In this appendix we collect various mathematical developments whose treatment in the body of 
the paper would interrupt the flow of the presentation. 

First of all we report several mathematical identities. We consider all of them well-known, 
but for completeness we either prove them below, or indicate where proofs can be found. These 
formulas feature sets of N numbers such as £„, or r] n or Cn.; these numbers are arbitrary but 
for simplicity we assume them to be distinct. The formulas of course remain valid when these 
numbers are not distinct, but possibly only by taking appropriate limits. Sometimes an arbitrary 
number z also appears. 



N 



E 

k=l 

N 

E 

k=l 

N 

E 

k=l 

N 

E 



N 



ii fci 



N 



Ci- 



k=l 



N 

n 



k n 

j=l, j^k 
N 

cr 1 n 

3=1, i+k 
N 

cr n 

— Cn 



1 



n r 



Cj - Cfc 



Cfc - Cj 



JnN, 



n = l,2,...,N, 



n=l,2,...,N, 



i=l, jfrn 



Cj Cn 



N 

E 

k=l 



N 

E 

k=l 



N 

n 



Cj ~ Vn 
0-6 



n, m = 1, 2 
Vj ~ Cfc 



? ■ ■ ■ 3 x ' 1 



N, 



N 

n 



N 



1 



z-Vk 



n (Cj-vk) 

j=l, ij£n 

N 

n (vj - vk) 

j=l, j=£k 



Vj - Vm 



N 



n, m = 1, 2, 



n (q-z) 

N 



N 

n 



Cj -z 



(A.l) 
(A.2) 
(A.3) 
(A.4) 
(A.5) 
N, (A.6) 



n = l,2,...,N, 



(A.7) 
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E 

k=l 



N 



n (d-m) 

j=l> j¥=n 



N 



n fa - vk) 

b=i, &k 



n = l,2,...,N, 



N 

E 

k=l 



N 

E 

fc=i 



1 



Vk ~ Cn 



Vk - Ck 

Vk- Z 



JV 



JV 

n 

3=1, jVfc 



JV 

n 



JV 

E 

fc=i 



JV 

E 

fc=i 



TV 

n (Cj-vk) 
j=i 

TV 

n - 

-i=i, jVfc 

TV 

% n (o - vk 

j=i, jV" 

TV 

n - %) 

- 3=1, j^k 



N 

E 

/<=! 



Vk - Ck 



JV 

n 



Vk - CnJ . " \Vj ~ Vk 

J=l, 3Tk 



Cj - Vk 



l 



^ - Vk 



Cj - Vk 

Vj - Vk 



N 

E 

k=l 



JV 



/<=i 



rc = 1,2,... ,A, 



(A.8) 
(A.9) 



JV 

n (Cj-vk) 

i=i 



/v 



(z - n (vj - v^ 

3=1, 3& 



(A.10) 



(A.ll) 



TV JV 
fe=l, fc^n fc=l 



n = l,...,iV. 



(A.12) 



The identity (jATTj) (with an arbitrary number) is implied by the fact that its left-hand side 
is a polynomial in z of degree less than N (in fact, of degree at most N — 1) which clearly has 
the value unity at the N points £ n , and the right-hand side, i.e. unity, is the unique polynomial 
of degree less than N in z that has the value unity in N distinct points. The identity (|A.2p 
is the special case of (|A.ip with z = 0. The identity (|A.3P coincides with equation (2.4.2-32) 
of [4] (or, as above, it is implied by the observation that its left-hand side is a polynomial in z 
of degree less than N the values of which at the N points Ck coincide with the values of the 
polynomial z n at z = Ck)- The identity ()A.4|) coincides with equations (2.4.3-12) and (2.4.3-21) 
of [I]. The identity (|A.5|) is obvious. The identities (|A.6p respectively (1A.7|) coincide with 
equations (2.4.2-26) respectively (2.4.2-27) of [I] (via the definition (2.4.2-24), with x n = £ n , 
Un = Vn, respectively x n = z, y n = r/ n , z n = Cn)- The identities (IA.8|) respectively (1A.9P follow 
from (|A.7p in the limit z — >■ oo respectively Cj ~^ °°- The identity ()A.10p follows from (|A.8P 
and (IA.7P via the trivial identity 



Cr. 



Vk 



Cn 



Vk 



Vk 



Finally, the identity (|A. 1 1[) follows from (|A. 10|) in the limit z — > oo, and the identity ()A.12p is 
just the special case of the preceding identity (jA.lip with Cn = 0. 

Next we report a simple lemma (for a neat proof see for instance [18], or below, after the 
proof of the following Lemma IA.2|) . 
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Lemma A.l. The solution of the set of N linear algebraic equations for the N variables gk 
reading 



N 

E 



9k 



is provided by the formula 

N 



n = l,...,iV 



(A. 13a) 



9k = c(£ k -r, k ) Y\ 

j=h i^k 



Cfe - Vj 



A generalization of this lemma reads as follows: 



fc = 1, . . . , JV. 



(A.13b) 



Lemma A. 2. The solution of the set of N linear 
reading 



N 

E 

k=l 



9k 



Cfe - Vn 



n = l,...,N 



rate equations for the N variables g k 



(A.14a) 



is provided by the formula 

N 



9k 



Vs 



=1 



N 

n 

J'=l, 3+k 



Zj - V* 

ii - ik 



N 

n 



n = l,...,N, 
or equivalently 

9k = (& ~ Vk) 
n = l,...,N. 



(A. 14b) 



N 

n 

3=1, j+k 



Vj - £k 



N 




r ^ 

n 


- 




E< 








j 




AT 






s=l 




n 

-i=i> j¥s 


- 





(A. 14c) 



To prove this formula one inserts this expression, ()A.14bp . of in (|A.14ap . and notes that 
one obtains thereby an equality provided there holds the formula 



N 

E 



€k - Vs 



or, equivalently, 

JV 



N 

n 



- % 



E 

fc=i 



N 

n 



N 

n 



N 

n 

i=i. j¥s 



^ - 6 

»7j - % 



A' 



Q ~ 7/Vi 



Clearly this formula is implied by the identity (|A.6p . Lemma lA.21 is thus proven. 

Note that, by setting c„ = c in (IA.14c|) and using the identity (|A,8|) (with the dummy index k 
replaced by s and the index n replaced by k), one reobtains (|A.13bp . thereby proving Lemma lA.ll 

We now report, and prove, 3 other lemmata. 
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Lemma A. 3. There holds the formula 

N 

det[I + X] = 1 + 

fc=i 



(A.15a) 



provided I is the N x N unit matrix and the N x N matrix X is defined componentwise as 
follows: 



X 



nun — *"m) 



n,m = 1,...,N. 



(A.15b) 



This (presumably well-known) formula is easily proven by recursion. 
Lemma A. 4. The N eigenvalues of the N x N matrix 

U nm — 0~nmCn ~i~ Imi 71, 771 — 1, . . . , iV, 

coincide with the N solutions of the following algebraic equation in z: 



(A.16a) 



N 

E 

k=l 



Vk 



z-Ck 



1 = 0, 



(A.16b) 



i.e. they are the N roots of the polynomial of degree N in z that obtains by multiplying the 

N 

left-hand side of this equation by Y\{z — Qk)- 

i=i 

This lemma is an immediate consequence of the preceding Lemma lA.31 because the secular 
equation associated with the matrix (|A.16a|) (whose roots provide the eigenvalues) is easily seen 
to coincide (up to an overall, hence irrelevant, multiplicative constant) with the vanishing of the 
determinant in the left-hand side of (|A.15ap with (|A.15b|) and Xf. = rjk/(z — Cfc)- 

Lemma A. 5. The inverse of the matrix defined componentwise as follows, 

fn9m 



Mr, 



n,m = 1, . . . , N, 



Cm Vn 

is defined componentwise as follows: 



(A.17a) 



[M- 1 ] 



9nfn 



N 

n 



N 

n 



Vj - in 



n, m = 1, . . . , N. 



(A.17b) 



The proof of this formula goes as follows. The matrix formula MM 1 = /, written compo- 
nentwise, reads, via ()A.17ap . 



N 

E 

fc=i 



'g k [M- r 



km 



Ck - Vn 



fn 



n,m = 1, . . . , N. 



Then, for fixed m, apply Lemma IA.2I with g^ replaced by gt[M 1 ] fcm , and c n replaced by 
Snm/fn- This yields, rather immediately, the formula (|A.17b|h which is thereby proven. 
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B Appendix 

In this appendix we detail the derivation of some findings for the fourth model, firstly the 
expression (|2.17bp of g n (z,z) as solution of the linear system (|3.35b|) . and secondly the derivation 
of (|2.19bj) from (|2.19a|) . 
Via the identity 



1 



1 



1 



V 



{j]z k + /3)(z k - az n - b) 7](az n + b) + (3 \z k -az n -b ijz k + P 

and the definition 

N 



a 



(0) 



E 

k=l 



V9k 



r\z k + PJ ' 



(B.l) 



(B.2) 



the linear system (|3.35bp can be conveniently reformulated as follows: 



N 

E 

fc=i 



9k 



z k -az n -b 



cn = a(0) + r 1 (az n + b) + p =a(Q) + a+ a 
r]z n + p Tj 



P{1 - a) 



az n + aP/rj 



We now use Lemma fA.21 (with £ k = z k , i] n = a z n + b and c n defined as above). We thus get 

N 



9n(z,z) = (z n - az n - b) 



n 

j=l, 3+ri 



z n — azj — b 

Z n Zj 



(1) ap{l - a) (2) 
V 



a 



(l) 



a 



(2) 



N 

E 

k=l 



N 

E 

k=l 



N 



(zj - az k - b) 

3=1, j^n 



N 



n = l,...,N, 



J] (azj - az k ) 
j=i, j= k 

N 

(zj - az k - b) 
1 J=l» r£n 

azk + ap/r] N 

[[ (azj-az k ) 

j=l, j=k 



n = l,...,N. 



It is now plain, via the identity (|A.8P (with 7] k = az k + b, Q = zj) that a„ = 1- As for 

(2) 

the sum er„ , it is also easily evaluated via the identity (IA.7P (now with r\ k = az k + b, (j = zj, 
z = b — aP/rj): 



a 



(2) 



-N„ N 



CL Tj 

rjz n + aP - br\ 1 \ \ r]Zj + P 

7 = 1 J 



n 



rjZj + a/3 — brj 



Hence 



9n(z,Z) = (z n - az n - b) 
N 

ff (3) = a i-^/3(i _ a ) Yl 



N 

n 

7]Zj + a/3 
i V VZj + P 



azj + b — z n 

Zj z n 



(3) 



■qz n + aP - brj 



(B.3a) 
(B.3b) 
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Now, using this expression of g n (z,z), we can evaluate do from its definition (|B.2p . thereby 
obtaining 



a 



(o) 



1 - a( 4 ) 



(3)^(5) 



-a + 



a + ^ 3 )q( 5 ) 
1 - 



a 



(l) 



a 



(5) 



N 

E 

fe=i 

N 

E 

fe=i 



z fc - az k -b 



N 



h + P/v 
Zk - azk-b 



n 



{riz k + P)(z k + a/3/r)) 



N 

n 



azj + b — Zk 

Zj z^ 



3=1, j^k 

It is now plain (via (|A.10p with ijj = Sj, Cj = a Zj + b, z = —fl/rj) that 
N 



aw=i-n 

3=1 



rjaZj + rjb + /3 
rjZj + /3 



(B.4) 



To evaluate 05, we use again an identity analogous to that used above: 



1 



1 



{j]z k + P)(z k + aP/rj) P(l - a) \z k + aP/rj r]z k + P 
Thereby 

0-(6) _ a (7) 



a 



(5) 



/?(1 - a) 
N 



a 



(6) 



tr 



(V) 



E 

fc=i 

N 

E 

fc=i 



4 -az k -b 



N 



n 

AT 

. n 



z-k + OLpjr] 
z k - az k -b 



azj + b — z k 

Zj z k 



azj + b — z k 

Zj z k 



(B.5a) 
(B.5b) 

(B.5c) 



Both these sums can be evaluated via the identity (|A.10|) . with r\ k = z k , Cj = azj + b, and with 
z = aP/r] respectively with z = —P/i], obtaining 



N 



3=1 

hence (via ()B.5aP ) 



rjZj + P 



rjZj + aP 



N 



CT 



(7) 



rjjazj + b) + p 
VZj + P 



a 



(5) 



1 



P(l -a) + b V 



N 

n 



Tjjazj + b) + p 
rjZj + P 



N 
3=1 



r/Zj + a/3 



and via this formula together with (|B.4p and (|B.3bp we finally get 



a^+a = a 



1N U 

3=1 



■qzj + a/3 
r/Zj + P 



The insertion of this expression of + a and of the expression (IB.3bj) of in (IB.3aj) 
completes the derivation of the expression (|2. 1 7b[) of g n (z, z) . 

Finally let us tersely outline the derivation of (|2.19bp from (|2.19ap . Firstly one uses in (|2. 19a|) 
the identity (jB.ip ; then one uses twice the identity (jA.lOp . with £ k = z k — b, r] k = az k and with 
z = a(az n + b) respectively with z = —aP/rj. And the rest is trivial algebra. 
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